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Abstract 

We carry out a detailed Lie point symmetry group classification of the Lienard type equation, 
x + f(x)x + g(x) = 0, where f(x) and g(x) are arbitrary smooth functions of x. We divide our 
analysis into two parts. In the present first part we isolate equations that admit lesser parameter 
Lie point symmetries, namely, one, two and three parameter symmetries, and in the second part 
we identify equations that admit maximal (eight) parameter Lie-point symmetries. In the former 
case the invariant equations form a family of integrable equations and in the latter case they form 
a class of linearizable equations (under point transformations). Further, we prove the integrability 
of all of the equations obtained in the present paper through equivalence transformations either by 
providing the general solution or by constructing time independent Hamiltonians. Several of these 
equations are being identified for the first time from the group theoretical analysis. 
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I. INTRODUCTION 



In this set of two papers we perform a Lie symmetry analysis for the Lienard type equation 

A(x, x, x) = x + f(x)x + g(x) = 0, (1) 

where over dot denotes differentiation with respect to time and f(x) and g(x) are arbitrary 
smooth functions of x. Notable equations from class ([T]) include a large number of physically 
important nonlinear oscillators such as the anharmonic oscillator, force-free Helmholtz oscil- 
lator, force-free Duffing and Duffing-van der Pol oscillators, modified Emden-type equation 
(MEE), and its hierarchy, and generalized Duffing-van der Pol oscillator equation hierarchy. 
These equations arise naturally in several physical applications. The outstanding represen- 
tative of the class of equations ([1]) is the modified Emden equation (also called Painleve-Ince 
equation), 'x + axx + /3x 3 = 0, where a and (3 are arbitrary parameters, which has received 
considerable attention from both mathematicians and physicists for more than a century 
(see for example Ref. [l| and references therein). 

During the past three decades, immense interest has been shown towards the search 
for symmetry generators of nonlinear ordinary differential equations (ODEs) and classifica- 
tion of low dimensional Lie algebras and linearization. Eventhough Lie himself had shown 
that the second order ODE of the form, x + f(t,x,x) = 0, can admit a maximum of 
eight symmetry generators, the recent impetus came only when Wulfman and Wybourne 2 
showed that the maximal Lie group of point transformations for the simple harmonic oscil- 
lator is eight and the associated group is SL(3,R). Subsequently, Cervero and Villarroel 3 
showed that the damped linear harmonic oscillator also admits eight symmetry generators. 
Thereafter, several studies were made to isolate the equations which admit rich Lie point 
symmetries by exploring their symmetry algebras and their applications in physics and 
mathematics^^^^^^^^. For a recent survey of results one may refer Ref. LLa and 
references therein. 

Specific equations of the form ([1]) have also been investigated from other points of view. 

rer symmetries for certain physically important systems were also studied 
Also contact symmetries for the harmonic oscillator were also explicitly 
constructed by Cervero and Villarroel^. The nonlocal symmetries for the MEE have also 
been studied in several papers^ 1 ^. Recently much interest has also been shown towards 



For example, Noet 

in Refs. HQQQ. 
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exploring generalized symmetries, namely, A-symmetries (also called c°°-symmetries) for 
certain nonlinear ODEs (see for example Ref. l2ll ). 

The main goal of the present set of papers is to present a detailed Lie point symme- 
try analysis of ([T]) and isolate integrable and linearizable cases explicitly At this point we 
mention here that the study of group classification is interesting not only from a purely 
mathematical point of view, but is also important for applications. Physical models are 
often constrained with apriori requirements to symmetry properties specified by physical 
laws, for example, from the Galilean or special relativity principles. In this work we not 
only isolate the invariant equations and symmetries but also present the integrals of motion 
and/or general solutions wherever possible. The motivation for the present study comes 
from our recent work in which we have considered a rather general second-order nonlinear 
ODE of the form x + (hx q + k 2 )x + k 3 x 2q+1 + k A x q+l + Ax = 0, where fcj's, i = 1,2,3,4, 
A and q are arbitrary parameters, and identified several interesting integrable cases and 
for certain equations we have also derived explicit solutions of both oscillatory and non- 
oscillatory types 22 . Interestingly, we demonstrated that a system which is very close to the 
MEE, namely 'x + axx + ^-x 3 + Xx = 0, possesses certain basic properties which is very 
uncommon to nonlinear oscillators^: It is a conserved Hamiltonian system (of nonstandard 
type) admitting amplitude independent harmonic oscillations. While the above mentioned 
equations are specific examples of (1), the question arises as to whether there exist other 
integrable/linearizable second order ODEs belonging to this class. In this and the accom- 
panying paper (referred to as II), we present a systematic analysis towards this goal. In this 
way we classify both integrable and linearizable equations which belong to the Lienard type 
system (1). 

Even though the Lie's algorithm is in principle a straightforward one, the group classi- 
fication for the present problem is reduced to integration of a complicated overdetermined 
system of partial differential equations for the infinitesimal symmetry functions for arbitrary 
forms of f(x) and g(x) in (1). While solving these determining equations we have to choose 
all the symmetry functions (see Eqs. ( |T2|) and ( fTBl) below) not equal to zero in order to 
obtain the maximal Lie point symmetries. On the other hand considering the special case 
of one or more of the symmetry functions to be equal to zero also, we obtain a spectrum 
of integrable equations. In this way we are able to classify (i) systems with non-maximal 
symmetries (less than eight) (ii) systems with maximal symmetries (eight). Note that for 
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second order ODEs the dimension of the Lie vector space cannot be four, five, six and seven. 
In the present paper we focus our attention only on the equations associated with lesser 
parameter Lie-point symmetry groups, that is one, two and three parameter Lie-point sym- 
metry groups. However, within this classification, we identify a wider class of important 
integrable equations which are invariant under 2-parameter point symmetry group. Many 
of them are being identified for the first time from the group theoretical analysis. In the 
second paper we identify all the equations admitting maximal number of symmetries which 
also turn out to be linearizable. The forms of / and g which lead to only two symmetry 
generators are as follows: 

(i) f = h + hx\ g = "j^-x + + (2a) 

X 2 fi A 2 A 2 

(ii) / = -fi + A 2 log(x), g = gix - + -j)x\ogx + -jx{\ogxf, (2b) 

C-\ t /l ^ , A h L 9 1 fey \ 

H/ = ^, 9 = ~T x + ^x- + -^ 2 (2c) 

(iv) / = £ + he-**, 9 = -^he~^ + 9l e-^ - |. (2d) 

(Here fx, g%, k 2 , Xi, A 2 , q, A are all constants). All the above equations are pointed out to 
be integrable through equivalence transformations. We also identify that the only system 
which admits a three parameter symmetry group is the Pinney-Ermakov equation where 
/ = and g = u?x — Jj (a;, g: constants). 

It is well known that the infinitesimal generators of a given Lie group form a Lie algebra. 
The Lie algebras constituted by Lie vector fields are widely used in the integration of dif- 
ferential equations^, group classification of ODEs and PDES— , in geometric control theory 
and in the theory of systems with superposition priniciples 2 ^ and in different schemes for 
numerical solution of differential equations 26 . A vast amount of works is available in the 
literature on the classification of realizations of finite dimensional Lie algebras on the real 
and complex planes. For example, the realizations of all possible complex Lie algebras of 
dimensions no greater than four were listed by Lie himself^. Recently Gonzalez-Lopez et al 
have provided the Lie's classification of realizations of complex Lie algebras 2 ^ and extended 
it to the real case. A complete set of inequivalent realizations of real Lie algebras of dimen- 
sion no greater than than four in the vector fields on a space of an arbitrary (finite) number 



of variables was constructed in Ref. 



29 



For more details on the classification of Lie vector 



fields one may refer the recent work of Ref. 



30 and references therein. On the other hand, in 



our paper we focuss our attention on constructing Lie vector fields for the class of equations 
resulting out of Eq. (1) alone and discuss their integrability. 

The plan of the paper is as follows. In the following section, we present the Lie's algorithm 
for Eq. (1) and discuss the solvability of the determining equations. A careful analysis of our 
investigations show that one should consider two separate cases, namely (i) the symmetry 
function (i) 6 = and (ii) b ^ 0, while solving the determining equations. Since the former 
case admits three symmetry functions, a(t, x), c(t, x), d(t, x), while classifying the integrable 
equations, we consider the possibilities (i) d = 0, a, c ^ 0, and (ii) c = 0, a, d ^ separately 
and bring out the equations that are invariant under both the possibilities in Sec. III. We 
also discuss sub-cases in both the cases (i) and (ii). Further, we prove that the system (1) 
does not admit a three parameter Lie point symmetry group when both f(x),g(x) ^ 0. 
In Sec. IV we investigate the equivalence transformations for Eq. (1) and show that they 
lead to integrable forms. In Sec. V we consider the special case in which either f(x) 
or g(x) is equal to zero and identify the associated integrable equations in this class. The 
notable example in this class includes Pinney-Ermakov equation. In Appendix A, we present 
some details on the Hamiltonian structure of an integrable equation that arises in the case 
(i). In Appendix B, we point out briefly some notable equations that are included in the 
most general equation (corresponding to (2a)). In Appendix C, we discuss the method of 
solving the integrable equation identified as integrable in this paper corresponding to (2b). 
The Liouville integrability of the two other integrable equations identified in the category 
c = 0, a, d 7^ are presented in Appendices D and E. Finally, we present our conclusions in 
Sec. VI. 

II. DETERMINING EQUATIONS FOR THE INFINITESIMAL SYMMETRIES 

We consider the one dimensional nonlinear Lienard type system of the form (CD). Let 
the evolution equation be invariant under the one parameter Lie group of infinitesimal 
transformations 

t = t + e£(t,x) + 0(e 2 ), 

i = x + er](t,x) + 0(e 2 ), e < 1, (3) 
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where £ and r] represents the infinitesimal symmetries associated with the variables t and x 
respectively. The associated infinitesimal generator can be written as 

X = + A. (4) 

Eq. (TjQ) is invariant under the action of (j4j) iff 

X^(A)\ A=0 = 0, (5) 

where 

is the second prolongation^ in which 

rjW = rj — x£, Tj = fj — x£ — 2x£, (7) 

and dot denotes total differentiation. By analysing Eq. (jSJ) we get the following determining 
equations: 

U = 0, (8) 
rh* ~ 26, + m* = 0, (9) 
2Vtx ~ £« + /& + 3g£ x + r)f x = 0, (10) 
Vtt ~ (Vx ~ 2&)</ + /fft + V9x = 0, (11) 

where subscripts denote partial derivatives. 
Solving Eqs. (jHJ) and we obtain 

£ = a (t) + b(t)x (12) 

and 

77 = te 2 -2^(s) + c(t)x + d(t), (13) 

where 

3 X = F(x) = / f(x')dx' and 9 M = f(x), (14) 
Jo 
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and a(t), b(t), c(t) and d(t) are arbitrary functions of t. With these forms of £ and 77, Eqs. 
(jTUj) and (jTTj) can be rewritten as 

(fee 2 - 2fe3 + cr + d)/ x + (a + bx)f + 3bg - AbF + 3bx + 2c - a = 0, (15) 

and 

(bx 2 - 2fe3 + cx + d)g x - (c - 2d - 2bF)g - 2b% 

+ (bx 2 - 2fe3 + cx + d)f + 'bx 2 + cx + d = 0, (16) 

respectively. Solving Eqs. (jT5j) and (jTBl) for the given forms of /(x) and g(x)we can get the 
infinitesimal symmetries. 

The foremost and simplest solution for ( JTBT) and (JTBT) for any form of / and ^ is a = 
constant, 6, c, <f = 0. In other words, one immediately gets the time translation generator 
X = irrespective of the form of / and g. Our motivation here is to find explicit forms of 
/ and g which admit more number of symmetries. For this purpose we solve Eqs. (1T51) and 
(flEj) in the following way. Rewriting Eq. f|T5l) . we get 

g = ^-{-(bx 2 - 263 + cx + d)^ - (d + fee)/ + AbF - 3bx - 2c + o], 6 ^ 0. (17) 

Thus the existence of Lie point symmetries of the general form (112j) and (|T3l) with 6^0 
introduces an interrelation between the functions / and g. However, this relation has to be 
compatible with the second determining equation f|T6|) . Thus using (ITTj) into (jTBT) . one can 
obtain an equation for /, (which also involves all the four symmetry functions) which fixes its 
form as well as the associated symmetries systematically. This is carried out in the following 
paper II, where we show explicitly that maximal (eight) number of Lie point symmetries 
exists only for the case f xx = and for f xx 7^ necessarily requires the symmetry function 
b = 0. Consequently, one has to consider the case b = in Eqs. (jT5j) - (jTBI) separately 
because of the condition b 7^ in Eq. (ITTj) . Thus it is of interest to consider two separate 
cases associated with Eqs. (jTBI) - (jTBI) : 

Case (i) b = : Since we assume one of the symmetry functions to be zero the determining 
equations lead us to lesser Lie point symmetries alone (one, two and three symmetries). 

Case (ii) b 7^ : In this case we solve the full determining equations which in turn lead us 
to the maximal (eight) Lie-point symmetry group (as well as other lesser point symmetries 
while b 7^ 0) when f xx = 0. 
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In this paper, we analyze in detail only the Case (i) and present the results of the other 
case in the subsequent paper II. 

A. Alternate Way 

One may also note here that one can proceed in an alternate way to analyze the deter- 
mining equations (|15|) and ([TBI) for compatibility to determine the forms of f(x) and g(x) 
and the associated symmetries. For example, Eq. (ITo*]) can be rewritten as 

a = (bx 2 - 263 + cx + d)f x + (d + bx)f + 3bg - AbF + 3bx + 2c. (18) 

Differentiating the above equation with respect to x and rearranging one can express b in 
terms of d, b, c and d. Differentiating again the resultant equation with respect to x and 
simplifying the latter we find 

df xx = -[(bx 2 - 263 + cx + d) xx f x - 2(bx 2 - 263 + cx + d) x f xx 

(bx 2 - 263 + cx + d)f xxx - bxf xx - 3bg xx + 2b f x . (19) 

Similarly from Eq. (Fl6|) one can obtain expressions for c and d and finally arrive at an 
expression of the form 

= -{(bx 2 - 263 + cx + d)g x - (c - 2d - 2bF)g - 263 

+ (bx 2 -2b^ + cx + d)f + "bx 2 ] xx . (20) 

To find a compatible solution of Eqs. (TT9"j) and (f2"0l one may substitute for d from (fT9"|) 
into (120]) and analyze the resultant equation to find the allowed forms of / and g and the 
associated symmetries. However, in practice we find the method leads to very lengthy and 
laborious calculations. On the other hand in the procedure we adopt we find that for the case 
fxx 7^ 0, one necessarily requires 6 = 0, see Sec. V in the following paper II. Consequently 
the analysis given in Sec. Ill follows naturally. 

On the other hand, for the case f xx = 0, Eq. (fl9|) leads to the condition 

= 2bff x -3bg xx . (21) 

which is consistent with (T2"U1) . From (|21[) one can immediately write f = fi + fi^ and 
9 — \ f1f2X 2 + |/| x 3 + 9i x + 92- To fix the corresponding symmetries one has to resubstitute 
the forms of / and g in the original determining equations and solve them consistently. This 
is carried out in Paper II. 
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III. LIE SYMMETRIES OF LIENARD TYPE SYSTEMS - LESSER PARAMETER 
SYMMETRIES: CASE b = 



We now explore the nature of the evolution equations which possess lesser parameter Lie 
point symmetries. Considering Eqs. (IT5|) and (TT5|) . we now assume the function b = to 
obtain the following determining equations, 



respectively. We note here that the determining equation (|22|) for the function / does not 
involve the function g. As a consequence an explicit form for / can be determined by 
direct integration. Now substituting this form of / into equation (123]) we can derive the 
corresponding form of g. 

It is a well known fact*^ that a second order ODE admits only 1, 2, 3 or 8 parameter Lie 
point symmetries (which we will see explicitly for Eq. (1) also in the present as well as in 
the follow up paper II). In Sec. 2 we noted that the most general equation which is invariant 
under the one parameter Lie point symmetry group is the general equation (pQ) itself with 
arbitrary form of f(x) and g(x) since there is no explicit appearance of t in the equation and 
the associated symmetry generator is But for b = 0, we explicitly show in the following 
that only two parameter symmetries exit when both / ^ and g ^ and obtain their 
specific forms, while three parameter symmetries can also exist only when / = 0, g ^ 0. 
Finally, in the case / ^ 0, g = 0, the second order ODE (1) can be rewritten as a first order 
equation which in turn can be integrated by quadratures straightforwardly. So we do not 
investigate this last category in this work. 

A. 2-parameter Lie point symmetries 




(22) 



and 



(cx + d)g x — (c — 2a) g + (cx + d)f + cx + d = 0, 



(23) 



Rewriting Eq. ( |22l) . we have 



Since / should be a function of x alone (vide Eq. (1)), we choose 



a , a — 2c , d , , 

- = Ai, = A 2 , - = A 3 , 25 

c c c 



where Ai, A 2 and A3 are constants. Then Eq. (1241) becomes 



(1 + A 3 ) (x + A3) 



Integrating Eq. fl26l) one obtains 

/ = ^ + /i(x + A 3 )- Al , (27) 

where fi is also an arbitrary constant. Note that since fi,X±, A 2 and A3 now occur in the 
form for f(x) which determine the ODE (1), they are not symmetry parameters but rather 
they are the system parameters. 

One can proceed further by constructing the associated form of g by solving Eq. fT23|) 
and classify the invariant equations. To start with, let us first classify for convenience the 
equations which are invariant under 2-parameter Lie-point symmetry group. To deduce 
these equations, we consider the two possibilities (i) d = 0; a, c 7^ 0, (ii) c = 0; a, d 7^ 
and also the sub-cases in both of them. We further note that due to the fact the system 
always admits translational symmetry, from Eqs. ([1]) and (j!2p . it is clear that a cannot be 
zero. So we need not consider the case a = 0, c, d 7^ 0. On the other hand if both c and d 
are simultaneously zero, then g = as may be inferred from Eq. ( |23l) . Finally, at the end of 
this section we consider the cases where none of the functions a, b, c are zero and show that 
even here only 2-parameter symmetries exist. 



1. Case 1 d = 0; a, c / (Ai + 0, A 2 + 0) 

The choice d(t) = with a(t), c(t) 7^ leads us to several interesting new integrable 
equations, as we see below. Solving Eq. ff25l) . with d(t) = and so A3 = 0, one can obtain 
explicit forms for the functions a and c as 



a = a 1 + ^(A 1 -2)c 1 e ( ^ )i , c = c 1 e ( ^ ) *, (28) 
A 2 

where a\ and c\ are two arbitrary (symmetry) parameters, which lead to a two parameter 
Lie-point symmetry group. Substituting Eqs. (I27j) and (128"]) into (123|) with d = 0, we get 

^ + x 9 ~ A 1 (A 1 -2)2 (Ai-2)* ■ [ ' 
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A 



Integrating Eq. (|29|) . we obtain 

P ~ A?(A 1 -2)2 X+ (2-A 1 )A 1 3; ' W 

where g\ is another integration constant. The above forms of / and g (vide Eqs. fl27|) and 
( 130|) . respectively) fix Eq. (pp) to the form 

^2 , ^-Ai^ , A 2( 1_ ^l)^, , Wl 1— Ai i 1— 2Ai 



* + u + /lX r + ik^w x + c*=w 1 + 911 = °- (31) 

For the sake of neatness, we rewrite Ai = —q, q 7^ 0, and A2 = —k 2 q, where k 2 is an arbitrary 
parameter, in the above equation so that we obtain 

x + (k* + A**) i + + J^f +1 + = °> ( 32 ) 

where k 2 , fi and q are nothing but system parameters. Eq. (1321) is the most general equation 
that is invariant under the two parameter Lie point symmetry group with infinitesimal 
symmetries 

fog hq 

£ = ai --(q + 2) Cl e(l + 2 ) . 77 = c x <M + 2 ) x . (33) 
k 2 



The corresponding infinitesimal generators read 

hq 



at 



+ x- 



&2 <9x 

The commutation relation between the vector fields X% and X 2 is given by 



(34) 



5. Integrability of Eq. (1321) for arbitrary values of g 



Eq. (1321) is the most general integrable equation which is invariant under the two pa- 
rameter symmetry group (I3"3"j) . Now we discuss the integrability of f[3"2"j) briefly here. By 
introducing the transformation 

w = z = e <9+ 2 > (36) 

qh 
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where w and z are new dependent and independent variables, respectively, one can transform 
(]32j) to the form 

w" + aw q w' + (3w 2q+1 = 0, (37) 

where a = J" 1 and 8 = ^jt 1 -. Eq. (1371) has been analyzed from different perspectives. 
For example, Lemmer and Leach 20 have studied the hidden symmetries of Eq. (1371) . Feix 
et al. 31 have shown that through a direct transformation to a third order equation the 
above Eq. (13"T|) can be integrated to obtain the general solution for the specific choice of the 
parameter 3, namely 3 = , . For this choice of 3, the general solution of (13"T|) can be 
written as 

x(t) = ( — — — -4— rf~~7 — n\~r~ I i h,h'- arbitrary constants. (38) 

For the same parametric choice recently we have shown that this equation can be linearized to 
a free particle equation through a generalized linearizing transformation so that the solution 
of the nonlinear equation can be constructed from the solution of the linearized equation^. 
However, our very recent studies show that Eq. (137]) admits time independent Hamiltonian 
description for all values of a and 3. By introducing appropriate canonical transformation 
to the Hamilton's canonical equation of motion one can integrate the resultant equations 
straightforwardly and obtain the general solution (for more details one may see Refs. 



34]). For convenience, the Hamiltonian structure of the above equation is indicated in the 
Appendix A. We also point out briefly other notable equations included in ( I32l) in Appendix 
B. 



3. Integrable equations with d = 0, a,c/0 (Ai = (q = 0), A2 ^ 0) 

Earlier while deriving the form (1271) for f(x) we assumed that Ai 7^ 0. Now let us consider 
the case Ai = 0. From Eq. ( 1251) we find that in this case 

-\ 2 t 

a — ax, c = C\e 2 , (39) 

where a± and C\ are two arbitrary symmetry parameters which again lead us to a two 
parameter symmetry group. Solving Eq. ( 124"1) . with the above forms of a and c, we obtain 

/(x) = -/ 1 + A 2 log(x), (40) 
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where fi is an integration constant. Substituting Eq. (140]) into Eq. (123]) with d = 0, we get 

<fe-§ + y(/i-A 2 logx) + ^ = 0. (41) 

Integrating Eq. (j4"Tl) . we obtain the following specific form for g, 

A 2 /i A 2 A 2 
9 = 9ix - (-|- + -j)xlogx + -jx(logx) 2 , (42) 

where #i is an integration constant. Using Eqs. ( 1401) and ( 1421) in Eq. (CQ), we have the 
following nonlinear ODE, 

A f A 2 A 2 
* + + A 2 log(a;))x + g x x - + -^)xlogz + -^x(logx) 2 = 0, (43) 

which is invariant under the following infinitesimal symmetries 

-A 2 t 

£ = ai, 77 = cie 2 x. (44) 

The associated symmetry generators take the form 

d >••>* d 

= . Txjj _. (45) 

The integrability of Eq. (|43p can be proved straightforwardly which we indicate in Appendix 
C. 



4- Integrable equations with d = 0, a, c ^ (Ai ^ 0, A 2 = 0) 

While deriving (T27|) we assumed that A 2 ^ 0. Now we analyse the case A 2 = with Ai ^ 0. 
In this case we find that the compatible solution exists for either Ai 7^ 2 or Ai = 2. In the 
first case by repeating the previous analysis we find that / = f\X~ Xl and g = ^x 3 ^ 1-2 j 
where fi and g\ are two arbitrary parameters so that Eq. (1) becomes 

x + f lX - Xl x + g!X {1 - 2Xl) = 0. (46) 

The associated infinitesimal generators turn out to be 

9 „ d x d . . 

Xl = Wf x > = % + (47) 



Eq. (146]) exactly coincides with (l3"Tj) by redefining Ai = — q, and so the integrability of (J46 
can be extracted from (13 
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In the second case, namely, Ai = 2, we obtain that the following form of equation for (1), 

fx. ^ Afi g\ 

T x + -7, ^ ~i 

x z 4 2x x A 



where A is an arbitrary parameter, which is invariant under the two parameter infinitesimal 
symmetry generators, 

*=£-*=«-"'(t£ + *I)- < 49 > 

We discuss the integrability of Eq. (l4"9l) in Appendix D. 

In this and previous sub-sections we discussed the cases (i) \\ = 0, A 2 7^ and (ii) A 2 = 
0, Ai 7^ 0. Finally, for the third case, namely (Hi) \\ = 0, A 2 = 0, one gets a = constant = 
di, c = constant = c\. The invariant equation turns out to be the linear damped harmonic 
oscillator equation 'x + j\x + g\X = 0, where fi and gi are arbitrary parameters. The 
associated infinitesimal vector fields are X\ = Jr, X 2 = xJj^. It is known that damped 
harmonic oscillator equation admits eight parameter symmetry group. Since one of the 
symmetry functions is zero we obtained only a two parameter symmetry group. The full 
symmetry group of the damped harmonic oscillator will be discussed in paper II. 

5. Case 2 c = 0, a, d / 0: Integrable equation 

In the previous sub-section we considered the case d = 0,a,c ^ 0. Now we focuss our 
attention on the case c = 0, a, d 7^ and fix the forms of / and g which are invariant under 
the corresponding symmetry transformations. Restricting to c = and a, d 7^ in (|22|) . we 
have 

+ 3' = \- (50) 

As before, since / has to be a function of x alone, we choose 



a a 

— = constant = Ai, — = A 2 = constant, (51) 
d d 



so that Eq. ( 150]) becomes 



Solving (|5T|) we obtain 



/x + A 2 / = A 1 . (52) 



An hl t hL t 

a = ai + —die*2 } d = die x ? , (53) 
Ai 
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where a± and d\ are two integration constants which are also the two symmetry parameters. 
Integration of Eq. fl52|) leads us to 



A 2 



A 2 x 



(54) 



where fi is an arbitrary constant. Substituting (1541) into (|23|) . with c = 0, we get 



A 2A 2 
g x + 2\ 2 g + -±f 1 e- x > x + -± = 0. 

A 2 An 



(55) 



Integrating fl55|) we obtain 



.9 



a 2 a 2 



(56) 



where gi is an integration constant. Eqs. ( 1541) and ( 1561) fix the equation (CTJ to the specific 
form 



x + 



A 2 



X X 2 
x — -rKjie + (?ie 



Af 



Af-- 



(57) 



We note that in the above Ai, A 2 , Qi, and fi are system parameters. Eq. (|57[) is invariant 
under the following two parameter Lie point symmetries 

An A l j 

£ = a(t) = a\ + —die^ , rj = d\e^ , (58) 
Ai 

where d\ and a\ are the symmetry parameters. The associated infinitesimal generators are 

d 



A| d_ d_ 
Ai <9i 9x 



(59) 



The integrability of the Eq. (1571) is discussed in Appendix E. 



B. Non-existence of 3-parameter symmetry group in the general case a, c, d ^ 

Finally, we consider the general case in which none of the functions a, c and d are zero. 
In this case, the function / takes the form given in Eq. (|27|) . The functions a, c and d can 
be fixed by solving the Eq. (1251) . Doing so we find 

Al(Al - 2) jj^t Tr^Sr* 7 \ TT^oT* /«n\ 

a = ai H cie( A i- 2 ) , c=cie< A i- 2 > , d = A 3 cie( A i- 2 ' , (60) 

A 2 

where only ai and C\ are the symmetry parameters. So in this case also only 2-parameter 
symmetries exist and no 3-parameter symmetry group is possible. 
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Substituting the forms /, a, c and d into Eq. (1231 and simplifying the resultant equation 
one obtains 

, (2Ai-l) _ 2A|(1-Ai) A 2 /i , . 

fc (s + A 3 ) * Ax(Ax - 2)2 + (2 - A0(x + A3)* ' 1 J 

where #i is an integration constant. One can directly integrate (|61j) to obtain 

a - AO (i + Xa) + ( V, A(» + a,)(;->.) + + Aa)(1 _ ai , 



" ' A,/ (2-A 1 )2 > " ' ' VAiy (2-A t 
Inserting the forms (1271) and (1521 in jT]) we get 

It is interesting to note that the system possesses only two parameter Lie point symme- 
tries. The infinitesimal symmetries and generators are 

g = 01 + Al ^ Al ~^ ciewV ^^'(i + As), (64) 
A2 

and 

9 v — i2_t Ai(Ai -2) d , XN <9\ 
Xl = -, X 2 = e^'(^4-J- + (x + A 3 y (65) 

respectively. 

Redefining X = X+A3 in (|63|) . the resultant equation coincides exactly with the integrable 
Eq. ( !3~TT) . The symmetry generators also coincide with the ones given in Eq. (j3~3~|) . So 
effectively no new nonlinear ODE is identified even when all the three symmetry functions 
are simultaneously nonzero. 

Thus we conclude that the system (1) does not admit a three parameter Lie-point sym- 
metry group when f(x),g(x) 7^ while the symmetry function b(x) = in (!H|)- (|TT|) . Further, 
the only equations which admit two parameter symmetry group alone are the four nonlinear 
ODEs given by Eqs. (13"2l) . (j43l . (j48|) and (ISTj). 



IV. EQUIVALENCE TRANSFORMATIONS 



We have shown in the above section that the identified evolution equations, namely Eqs. 
( 132|) . (}4"3"j) . (j48|) and (1571) . admitting two parameter Lie point symmetries can be transformed 
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into integrable equations (1321) . (IClj) . (IDlj) and (lElj) respectively, through appropriate trans- 
formations. In this section we give a group theoretical interpretation for these results through 
equivalence transformations (ETs). We invoke the equivalence transformations and give an 
explanation for the results since the group classification problem is closely related to the 
concept of equivalence of equations of the above forms with respect to transformations, see 
for example Ref. [32. 



Considering our original differential equation (1), let us consider a set of smooth, locally 
one-to-one transformations T : (t,x,f,g) — > (T,X, f%, g\) of the space M 4 that act by the 
formulas 

T = F(t,x), X = G{t,x), f 1 = H(t,x,f), 9l = L{t,x,g) (66) 

A transformation is called an Equivalence Transformation (ET) of the equality x = —f(x)x — 
g(x) if it transforms the equation 

x = —f(x)x — g(x) (67) 

to an equation of the same form 

X = -f 1 (X)X-g 1 (X). (68) 

In this case, Eqs. ( IBTj) and ( 1681) and the functions {f(x),g(x)} and {f\(X),gi(X)} are 
equivalent^. 

It is a proven fact that equivalent equations admit similar groups (for local transforma- 
tions) and ET is a similarity transformation. That is, if (IST)) admits the group E then (|68|) 
also admits a group similar to it for local transformations. 

Substituting the transformation (I6"6"|) into Eq. (I6"8"|) we get 

{G 2 t F t h + G 3 t9l ) + x[(G 2 t F x + 2G t F t G x )h + 3G 2 t G x9l ] + x 2 [3G 2 x G t9l 
+ {F t G 2 x + 2G t F x G x )h] + x 3 [F x G 2 x h + G 3 x9l ] = -{G t + xG x )[{F u 
+ 2xF tx + x 2 F xx - F x {xf + g)} + (F t + xF x )[(G a + 2xG tx 
+ x 2 G xx -G x (xf + g)l (69) 

where the subscripts denote partial derivative with respect to that variable. Equating the 
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coefficients of different powers of x n , n = 0, 1, 2, 3, we get 

FxGxfi + G x g\ = F X G XX — G X F XX , (70) 

(F t G 2 x + 2G t F x G x )f 1 ] + 3G x G t gi = F t G xx + 2F x G tx — G t F xx — 2G x F tx , (71) 
{G 2 F X + 2G t F t G x )f 1 + ?>G 2 t G x g\ = —2G t F tx — G x F tt + fF x G t + 2F t G tx 

- F t G x - F x G tt , (72) 

G 2 F t fx + G 3 t9l = -G t F tt + gF x G t + F t G tt - gG x F t . (73) 

Solving Eqs. ( 1701) and ( 1TT1) consistently we find G x = and F xx = 0. As a result one gets 

G = a(t), F — (3(t)x + 7(t), (74) 

where a, /3 and 7 are arbitrary functions of t. Substituting Eq. (ITlj) in (IT2"j) and (1731 and 
simplifying the resultant equations we get 

a 2 fij x = f{3a + P&- 2a/3, (75) 
a 2 ($x + 7)/x + d 3 #i = — a{j3x + 7) + g /3d + d(/^ x + l)- (76) 

From Eq. fl7"5l) we can obtain an expression which connects the transformed function f\ with 
the original function / of the form 

a a z pa 

Substituting (1771) in fl7"6l) and simplifying the resultant equation we arrive at 



Pg (Px+j) f , 2/j(/3x + 7 ) (/3x + 7 ) , 7S , 
cr cr per cr 



Thus we obtain the general ET 



T = a{t), X = /3(t)x + >y{t), f 1 = L + ?L-W 

a cr pa 



ft? (^x + t) 2^x + t) (/fo+j) ^ 

3i = 9i = To — 2 / H tttq • I '9) 

cr cr per cr 



Since we have already identified only four equations (vide Eqs. (132]) . f|43|) . PB]) and flSTl) ) 
that are invariant under two parameter Lie point symmetries within the class of equations 
(1) we consider only these four equations and present our result. Now solving Eqs. ( 1791) 
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with the given form of / and g one obtains the following result 

( Q -L 2) 1 k 2 k 2 

Case 1 (Eq.([32])) a = - , ; e Wt, (3 = e Wt, 7 = 

so that h = aX q g x = (3X 2q+1 
Case 2 (Eq.(g3D) a = t, (3 = e -^/M*M]* 7 = 

so that fi = constant gi = X 
Case 3 (Eq. (ggj)) a = -e A \ (3 = e At , 7 = 

so that /i = — 0i = 4g 

Case 4 (Eq.([57D) a = -^e"^, p = 1, 7 = -^ 

Ai A 2 

so that A = e A2f/ ^ = e" 2A2C/ (80) 



It directly follows that with the above form of f\ and gi, Eq. ( 1681) takes the form of ( |32l ). 
(ICip . (ID lh and (lEip ) respectively, which were shown to be integrable. 

V. LIE SYMMETRIES OF EQ. (pQ) WITH f(x) = OR g{x) = 

Next we consider the special case of Eq. ^ with f(x) = 0, that is, 

x + g(x) = 0. (81) 

In the following we focuss our attention only on the case b = so that we have £ = 
a(t), i] = c{t)x + d(t). Eqs. (TlOT) and (TTTT) with b(x) = and /(x) = give rise to the 
following conditions, respectively, 

a - 2c = 0, (82) 

and 

/2a — c\cx + d . , 

^+7^0+7^ = °- (83) 



-ex + d/ cx + d 
Since ^(x) should be a function of x alone, we choose 

--l = A l5 - = A 2 , ^ = -A 3 , (84) 
c c c 

d c 

where Ai, A2 and A3 are constant parameters. Note that the above implies - = A2- = A2A3. 

c c 
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Solving (154")) we find that the solution exists either for the parametric choice Ai 7^ 3 or 
Ai = 3. The respective infinitesimal symmetries are 

«c = a (t) = at + a 2 t, rj = c(t)x + d(t) = — + - Ai ^ 3 (85) 

[1 + Ai) [i + Ai) 

£ = a {t) = ai - - t + ^l\ c 2 cos y%t - ci sin a/A^*), 
2v A3 



?7 = c(t)x + d(t) = a± — (ci cos y X^t + c 2 sin a/ Ast)(a; + A2), Ai = 3 (86) 
The respective invariant equations turn out to be 

i+ (?TkF = ' Al ^ 3, (87) 

and 

x + ^(x + A 2 )+ l\ = 0, A 1 = 3. (88) 
4 (2 + A 2 ) ,:S 

Thus Eq. (|87p admits a two parameter symmetry group with the generators 



X - 9 X -t 9 I 2(X + A2) ^ (891 

On the other hand Eq. ( 1551) admits a three parameter symmetry group with the symmetry 
generators 



^»A'llf5 + (l + W il' (90) 



Redefining rc + A 2 = X in Eqs. (1571) and (1551) we get 

X + ^ = 0, A^3. (91) 

X + u; 2 X-^ = 0, ^ 2 = ^, £ = -<7i, Ai = 3, (92) 

Equation in (|91~T) corresponds to a conservative Hamiltonian system (if = y + ^ 1 ^^ X 1 ~ Al ) 
and so the Liouville integrability is assured. On the other hand Eq. (1921) is nothing but the 
Pinney-Ermakov equation, whose origin, properties and the method of finding its general 
solution have been discussed widely in the contemporary nonlinear dynamics literature (see 



for example Ref . 



33l and references therein) . For the sake of completeness we give the general 



solution of this equation as 



X = ~^\/(^ ~ 9) sm 2 M + <f>) + 9, (93) 
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It has also been shown that Eq. (1921) can be transformed to harmonic oscillator equation 
through suitable nonlocal transformation and from the solution of the latter one can con- 
struct the solution for the nonlinear equation. For more details one may refer—. 
Finally, for g(x) = 0, Eq. ([T]) can be written as 

x + f{x)x = 0. (94) 

Eq. (1M1) can be transformed to a first order equation by a trivial change of variable which 
in turn can also be integrated trivially. So we do not discuss the symmetries of this equation 
here. 

VI. CONCLUSIONS 

In the present paper we have investigated the Lienard type equation (1) in the framework 
of modern group analysis of differential equations. Even though the integrability properties 
of some of the specific equations coming under the Lienard type have been discussed in the 
literature, we have identified all those equations which admit only two and three parameter 
symmetry groups. 

To identify the integrable equations belonging to the class (1) we have deduced all the 
equations that are invariant under one, two and three parameter Lie point symmetries. Ob- 
viously the general Eq. does not contain the variable t explicitly and so it always admits 
a time translational generator. However, we have demonstrated that several equations admit 
two parameter Lie point symmetry groups. In particular these equations correspond to four 
specific forms of the functions f(x) and g(x) in (1), see Eq. (2), namely Eqs. (1521 . (1351) . 
( HHj) and (1571) . These equations have been deduced here through a group theoretical point 
of view alone. We have also discussed the integrability properties of these equations briefly 
and shown the existence of equivalence transformations. After analyzing the Lie point sym- 
metries we have also shown that Lienard type equation does not admit a three parameter 
symmetry group when both f(x),g(x) ^ in Eq. (1). However, in the sub-case, f(x) = 0, 
one can find that the well known Pinney-Ermakov equation is the only equation which is 
invariant under a three parameter Lie point symmetry group. 

In this paper we have restricted our attention only on the non-maximal Lie point symme- 
try groups. The question which naturally arises is what happens if one considers the more 
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general case, b ^ 0, vide Eq. (ITT}) . Such an analysis allows us to isolate a class of equations 
admitting eight parameter symmetries. We will present the results in the follow-up paper 
II. 
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In the following, we briefly discuss the integrability properties of the equations derived 
in Sec. III. To begin with let us consider the Liouville integrability of Eqs. (33). 



APPENDIX A: TIME INDEPENDENT HAMILTONIAN FOR fl37j) 



Recently, we have studied the integrability of (132]) or equivalently (157)) and found that it 
admits time independent integrals for all values of the parameters a and / 3 22 ' 34 . From the 
time independent integrals we have identified the following time independent Hamiltonian 
for (157)) . namely, 



H 



(r-l) !r=*l 
(r-2)^ 



(f— 1) * o+l 

- — '-apw 1 ^ , 



\pw q+l + log( 



a 2 > A(3{q + 1) 
a 2 = A(3{q + 1) 



- %pw q+l , a 2 < 4/3(g+l) 



where the corresponding canonically conjugate momentum is defined by 

(l-r) 

j [w + ^aw" +L \ 

P 



(9+1) 



„1+ 



r tan 1 



awi +1 +2(q+l)w 
2UJW1+ 1 



a 2 > A/3{q + 1) 
a 2 < 4/3( g +l), 



where r = 2 p(q+i) ( a ^ \A* 2 ~ 4/3(g + 1), u = \\J Aj3(q + 1) — a 2 and a = - 
details about the derivation of the above Hamiltonian one may refer to Ref. 
independent Hamiltonian ensures the Liouville integrability of (152)) or (1571) . 



(Al) 



i 



(A2) 



For more 



34. The time 
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APPENDIX B: NOTABLE INTEGRABLE EQUATIONS IN (|521) 



Besides the general case, q = arbitrary, Eq. (1321) encompasses several known integrable 
equations of contemporary interest. The interesting equations can be identified by appro- 
priately choosing the parameter q as we demonstrate briefly in the following. 

For example, choosing q = 1 in (132]) one gets the generalized MEE, 

2k 2 k f 

x + (k 2 + f x x)x + — ^x + -^V + gi x 3 = 0, (Bl) 
9 3 

Eq. ( 1B1I) can be transformed into the MEE, w" + f\ww' + giw 3 = 0, by introducing a 
transformation iy = xe~3~ and z = —-^-e~s~ . The Hamiltonian structure for this equation 
can be extracted from (lAlj) by restricting g = 1 in the latter relations. The restriction 
fi = in ( 1B11) provides us the force-free Duffing oscillator whose invariance and integrability 



properties have been discussed in Refs Jl4jj35l . With the choice f\ — 3, g\ — \,k 2 = 0, the 



resultant equation becomes a linearizable one whose invariance and integrability properties 



36 



37 



38 



39. 



have been discussed in detail in Refs. 
The case q = 2 in (1321) gives us 

3k 2 k f 

x + (k 2 + f lX 2 )x + — ^x + ^fV + gi x 5 = 0. (B2) 
lo 4 

The explicit form of the Hamiltonian can be fixed from (lAlj) by restricting q = 2 in the latter 
relations. We note here that Eq. (1B2|) also includes several known integrable equations. 
The notable examples are force-free Duffing-van der Pol oscillator equation (gi = 0) and the 
second equation in the MEE hierarchy (/i = 0,gi — ^). 

Finally, we note that one may also recover specific equations like the force-free Helmholtz 
oscillator and the associated Lie symmetries can be obtained by appropriately choosing the 
value of the parameter q. Choosing q — | and fi = in (1321) one gets the force-free Helmhotz 
oscillator. The symmetries of (1331) with q = | coincide exactly with the one reported in Ref. 



40. 



APPENDIX C: METHOD OF INTEGRATING EQ. d43]) 

The solution of Eq. (l4"3l) can be constructed from the solution of the damped harmonic 
oscillator using a general procedure given by us sometime ago in Ref. I4I. For example, let 
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us consider a linear ODE of the form 

U + aU + gi U = 0, (CI) 

where a and g\ are arbitrary parameters. By introducing a nonlocal transformation of the 
form U = xe~ 

■ J io g z(f)df in the Hnear ODE (JUT]) the latter can be brought to the form 

A ex A A^ 

x + (a + — + X 2 \og(x))x + g\x + log x + -^x(log x) 2 = 0. (C2) 

Now redefining the constants (at + ^) — ~.fi in flC2p one exactly ends up with fT4"5j) . 

Following the procedure given in Ref. Id one can obtain the general solution for ( 1C2I) 
from the linear equation. 



APPENDIX D: METHOD OF INTEGRATING EQ. (|48j) 



Eq. (jUJ) can be transformed to the equation of the form 

0. 



u 2 u 3 



(Dl) 



through the transformation U = xea*, Z = \e At . Eq. (1D1I) admits Hamiltonian structure 
for all values of fi and g\. The underlying Hamiltonian reads 



H 



(r-l). 



(r~2) 



I/' 



i<n(r-l) _|_ 

^ + §log [^sec^], 
where the canonical conjugate momentum is defined by 



P 



_ (r-l)/n (l-r) 



=2 tan" 1 [^2] 



fi > 

/i 2 = -4<7i 
A 2 < "40!, 



A 2 > 4/7! 
/i 2 < 4<7i 



;D2) 



(D3) 



where r = -^(/i ± V/ 2 + 4^) and w = ± ^/-4 5l - / 2 . 

The time independent Hamiltonian given above ensures the Liouville integrability of Eq. 

dm}. 



APPENDIX E: METHOD OF INTEGRATING EQ. ([5711 

By introducing a transformation U = x — and 2 = -^-e - ^* in (1571) the latter can be 
transformed into the form 

U" + he^ u U' + 3ie- 2A2t/ = 0. ('= ^-) (El) 
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Eq. (IEip can be rewritten in the form form 



U" + hf(U)U' + gJ(U) J f(U)dU = 0, (E2) 

where f(U) = e~ X2U and g± = — \2gi- Eq. flE2j) admits time independent Hamiltonian for 
all values of f± and g±. The respective Hamiltonians are 



H 



2A 2 



log[p] + irr-c 



|l0g 



A, 



sec 



-A 2 



2A 2 ^ C ' 



/i 2 > 4ft 
fx = *9i 
fi < 4^i, 



(E3) 



where the canonically conjugate momentum is defined by 



P 



[U> 



r A 2 J 



-A 2 e A 2^ t -1 r 2A 2 t/'-/ie- A 2^ 1 



2we 



/i 2 > 4^i 



(E4) 



where r = ^(^ ± ^ -4ft), w = ±x/4ft - Jf. 

One may note that the above Hamiltonian resembles the Hamiltonian structure of (I AID . 
The reason for this is that both the Hamiltonians flE3p and (1 A 1 [) can be generated from 
the time independent Hamiltonian of the damped harmonic oscillator by suitable nonlocal 
transformation. For more details about this nonlocal transformation one may refer to Ref. 
34. 
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